Sensors 2009, 9, 7988-8006; doi:10.3390/s91007988 



OPEN ACCESS 



sensors 

ISSN 1424-8220 

www.mdpi.com/joumal/sensors 

Article 

A Coupled Field Multiphysics Modeling Approach to 
Investigate RF MEMS Switch Failure Modes under Various 
Operational Conditions 

Khaled Sadek, Jonathan Lueke and Walied Moussa * 

Department of Mechanical Engineering, University of Alberta / University of Alberta, Edmonton, AB, 
T6G 2G8, Canada; E-Mails: kmahmoud@ualberta.ca (K.S.); lueke@ualberta.ca (J.L.) 

* Author to whom correspondence should be addressed; E-Mail: Walied.Moussa@ualberta.ca; 
Tel.: +1-780-492-6027. 

Received: 30 July 2009; in revised form: 15 September 2009 / Accepted: 16 September 2009 / 
Published: 12 October 2009 



Abstract: In this paper, the reliability of capacitive shunt RF MEMS switches have been 
investigated using three dimensional (3D) coupled multiphysics finite element (FE) 
analysis. The coupled field analysis involved three consecutive multiphysics interactions. 
The first interaction is characterized as a two-way sequential electromagnetic 
(EM)-thermal field coupling. The second interaction represented a one-way sequential 
thermal-structural field coupling. The third interaction portrayed a two-way sequential 
structural-electrostatic field coupling. An automated substructuring algorithm was utilized 
to reduce the computational cost of the complicated coupled multiphysics FE analysis. The 
results of the substructured FE model with coupled field analysis is shown to be in good 
agreement with the outcome of previously published experimental and numerical studies. 
The current numerical results indicate that the pull-in voltage and the buckling temperature 
of the RF switch are functions of the microfabrication residual stress state, the switch 
operational frequency and the surrounding packaging temperature. Furthermore, the 
current results point out that by introducing proper mechanical approaches such as 
corrugated switches and through-holes in the switch membrane, it is possible to achieve 
reliable pull-in voltages, at various operating temperatures. The performed analysis also 
shows that by controlling the mean and gradient residual stresses, generated during 
microfabrication, in conjunction with the proposed mechanical approaches, the power 
handling capability of RF MEMS switches can be increased, at a wide range of operational 
frequencies. These design features of RF MEMS switches are of particular importance in 



Sensors 2009, 9 



7989 



applications where a high RF power (frequencies above 10 GHz) and large temperature 
variations are expected, such as in satellites and airplane condition monitoring. 

Keywords: RF MEMS switch; substructuring; skin effect; buckling; residual stresses 



1. Introduction 

Radio Frequency (RF) Micro-Electro-Mechanical Systems (MEMS) capacitive shunt switches have 
had an enormous impact on the applications of wireless communication systems and satellite 
technology. These switches raise the performance of communication systems to a level significantly 
above their solid-state counterparts such as FETs and PIN diodes. The prime reasons that these RF 
MEMS switches are currently under immense consideration includes, their wider capacity frequency 
range (from RF to mm-wave signal), lower power consumption, lower insertion loss, and lack of 
inter-modulation distortion [1-4]. The reliability of these switches are reported to sustain up to 0.1^0 
billion cycles, while only handling power with maximum capacity of 0.5 W [4-7]. The major failure 
modes of RF MEMS switches are shown to occur due to membrane stiction, buckling, and creep [5]. 

The operation of RF MEMS switches requires an electrostatic force in order to actuate from the up, 
or "off position, into the down, or "on" position. The voltage required to actuate the RF MEMS 
switch in this manner is referred to as the pull-in voltage. This voltage is dictated by the geometry and 
stiffness of the RF MEMS switch. A higher pull-in voltage will decrease the overall lifespan of the RF 
MEMS switch, as shown by [1]. 

As in all micromachined MEMS devices, residual stresses develop in the manufacturing process. 
The residual stresses can be separated into two specific stresses, thermal and intrinsic stresses [8-10]. 
Thermal stresses are generated through differences in thermal expansion coefficients in adjacent 
materials in the manufacturing process. The generation of intrinsic stresses is far more complicated. 
Generally, intrinsic stresses involve any stress generated from excess vacancies, crystal dislocations, 
grain boundary interactions, or phase transformations [8,10]. Overall, the residual stresses generated in 
RF MEMS switches ultimately affect the general state of stress, and therefore, failure mechanisms of 
the MEMS device. While loaded with high frequency alternating current, the RF MEMS switch 
undergoes an electron crowding phenomena, otherwise known as the skin effect [5]. Due to self- 
inductance, the electrons flowing through the switch will be crowded on the edges of the conductor, 
creating a localized area of high temperature near the edges of the RF switch. Therefore, severe 
compressive thermal stresses develop due to this localized heating, causing the switch to buckle 
upwards. Depending upon the magnitude and direction of the residual stresses, this localized heating 
could potentially cause the RF MEMS switch to fail under severe buckling with sufficient load. 

The geometry of the RF MEMS switch greatly affects the power handling and frequency loads that 
each switch can experience. In addition, the actuation voltage of RF MEMS switches is affected by the 
geometry of these switches; for example, a mechanically stiffer switch will require a larger pull-in 
voltage. Mechanical approaches such as the addition of corrugations and holes can be introduced to 
control the membrane stiffness of the RF MEMS switch, in order to attain a specific design pull-in 
voltage [4]. 
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In this paper, a multi-domain 3D FE model is developed to investigate the multiphysics interactions 
within capacitive shunt RF MEMS switches. This coupled field model studies the electromagnetic, 
thermal, structural and electrostatic physics of RF MEMS switches, in response to variations within 
geometry, residual stresses (mean and gradient) and operational current frequency. The obtained 
results are then utilized to predict the onset of stiction and buckling with RF MEMS switches and its 
effect on these switches long term reliability. 

2. Coupled Field FE Analysis of RF MEMS Switch 

A flowchart illustrating the analysis procedure of the RF MEMS switch coupled field FE model, is 
shown in Figure 1. The full model was solved numerically using ANSYS 10.1 indirect multiphysics 
solver [11]. 



Figure 1. Flow chart for coupled field analysis procedure for RF MEMS switch. 
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For the current study, the RF switch model consists of corrugated aluminum membrane, silicon 
substrate and bottom electrode, as shown in Figure 2. 

First, the distribution of the current density (J) is calculated at a frequency (od) range of 0.1-100 GHz 
by conducting a full wave harmonic electromagnetic (EM) analysis. The heat generated in the switch 
membrane is calculated from J, and then used to calculate the temperature distribution (7) through a 
steady state thermal analysis. The EM-thermal analysis was performed as a sequentially coupled field 
analysis involving two-way multiphysics coupling (i.e., iterations are performed between both fields). 
Second, a coupled thermal-structural analysis is performed to calculate the in-plane thermally induced 
stresses and resulting deformation in the switch geometry. The coupled thermal-structural analysis 
involves one-way multiphysics coupling between both fields. Finally, the pull-in voltage {Vpi) for the 
new switch configuration, was calculated using a two-way coupled structural-electrostatic FE analysis. 
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Figure 2. Corrugated RF MEMS switch. 
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2.1. EM-Thermal Fields Coupling 

In the full wave EM field, a current (7) load that varies harmonically with time, is applied at the 
edge (x=0) of the switch membrane (Equationl). The electric field distribution is determined by 
solving the governing Maxwell's equations with two known BCs (Equation 2): 

I = I^^m{(Dt), x = 0 (1) 



dB dE 
V E = , VB = 0, and B= u J + s u — 



BCs, E = \^^ , B.s = 0 



(2) 



where, lo is the current amplitude; V, is the divergence operator; E, is the electrical field strength; 
Pm(T), is the temperature dependant electrical resistivity of the switch membrane material; So, is the 
permittivity of free space; V^, is the curl operator; 5, is the magnetic field strength; and //^ is the 
permeability of free space. 

After solving for E, the current density is calculated by Equation 3, as follows: 

p7^) 



Once the current density is calculated from the EM field, the model is switched to the thermal field. 
The heat generation {H) throughout the beam is calculated by Equation 4: 



H = \\J\\ PJT) 



(4) 



For the frequencies considered in this study (co > 0.1 GHz), the time constant for the thermal 
response is much longer than the period of variation of the input power. Therefore, the temperature 
distribution throughout the switch membrane is calculated using the steady state heat equation given 
by [5]: 
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VikVT) = -\\jf p^iT) (5) 

where, km is the thermal conductivity of the switch membrane material. The convection and radiation 
heat transfer mechanisms were negligible compared to the heat transfer by conduction. Therefore, only 
conductive heat transfer was considered in the thermal analysis. Same assumption has been validated 
in a previous study reported by Jensen et al [5]. 

For the thermal model, the boundary conditions (BCs) assume that the edges of the switch 
membrane (x = 0, x = Im) are kept at room, or in other words the reference temperature (Tj), Adiabatic 
BCs are assumed for the other four sides of the membrane (z = go. z = go + tm, y = 0, midy = Wm)^ 

Once the temperature distribution is found, a new value for the electrical resistivity is calculated and 
the analysis is switched back to the EM field to update the current density distribution. Iterations 
between the EM and thermal fields are executed in sequence until a convergence for the value of T is 
achieved. The current applied to the EM field is ramped up and the whole procedure is repeated until 
the value of the maximum current is reached. 



2.2. Thermal-Structural Fields Coupling 

The temperature distributions calculated from the coupled EM-Thermal FE analysis are applied to 
the structural field model as body loads. The coupling between the thermal and structural fields was 
conducted in a one-way sequential interaction manner. The model was completely switched to the 
structural field after the thermal analysis has been fully conducted. No iterations were performed 
between the two fields. Due to the mismatch of the coefficient of thermal expansion between the 
switch membrane and the substrate, thermal stresses and a new switch deformation state results from 
the applied T. The thermally induced axial strain (8th,x) and stress (ath,x) are calculated at a given 
position (x, z) by [5]: 

T{x,z) 



where, Aa is the difference in thermal coefficient expansion between the membrane and the substrate; 
Em is the effective Young's modulus of the switch membrane, given by Em = Em/il-Vm^) for 
Wm>5tm[l2]; Em and v^, refer to the membrane material Young's modulus and Poisson's 
ratio, respectively. 

Once the thermal induced stress distribution is evaluated, the new switch geometry is found by 
calculating the membrane deflection (Dm) from Karaman equation, which can be expressed in 1-D 
form as follows [13]: 



d'D 1 

m_ 



E t fdD ^V^2n A 



dx 



d'D . 
dx^ 



n 

dx^ 



=0 



and — ^ = 0 

dx i^-u.-Lffi 



(7) 



where, D is the flexural rigidity of the membrane, given by: 
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(8) 



The first term in right hand side of Equation 7 represents strain stiffening, and is considered only 
for large membrane deflections [12,14]. 

2.3. Structural-Electrostatic Fields Coupling 

The non-uniform thermal induced stresses calculated from thermal- structural coupled field analysis 
results in a new configuration for the switch geometry, as shown in Figure 3. The new geometry is 
used as a starting point for the structural-electrostatic analysis. Initially, mean and gradient residual 
stresses measured during microfabrication [8], were applied to the switch to account for the final 
geometry and state of stress before proceeding to the coupled structural-electrostatic analysis. 

Figure 3. Distribution of displacement vertical component along x-coordinate predicted by 
coupled EM-thermal analysis followed by a coupled thermal-structural analysis for 
corrugated aluminum RF MEMS switch membrane (/^ = 91 jum, = 20jum, t^ = OAjum, 
tc = 4tm, Nc = 13, CO = 20 GHz, P = 0.75 W). 




In the electrostatic field, the essential BCs are the applied ramped voltages at the bottom and top 
conductors. Since the width (Wm) of the modeled switch membrane is greater than the initial gap length 
(go), the fringing field effects were neglected [5,12,15,16]. The voltage distribution was calculated by 
solving Laplace equation with two known BCs [13]: 



d^V d^V d^V 



= 0, 



dx' dy' dz' 



(9) 



and F = 0 
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where V refers to the voltage distribution and Vo is the appUed ramped voltage at the bottom of the 
aluminum switch membrane. 

The electrostatic force distribution was then calculated, and the model is switched to the structural 
field. A new deformation state was obtained from the structural analysis and the model was switched 
back to the electrostatic field. Morphing BCs were applied at the side and top boundaries of the 
underlying medium to account for the new deformation state of the switch [11]. The electrostatic field 
mesh was updated and a new value for the electrostatic force distribution was calculated. Iterations 
between the structural and electrostatic fields are executed in sequence until a convergence was 
reached. The applied voltage is ramped up and the process is repeated for the new voltage value. The 
pull-in voltage (Fp/) was finally identified when a convergent solution was not achievable [12,13]. 

3. Automated-Substructuring Algorithm 

An automated substructuring algorithm was developed, in the current study, to reduce the 
computational cost of the coupled multi-field RF switch model. The same algorithm was utilized, in a 
previous study reported by the authors, to solve the coupled electro-thermal FE problem of MEMS gas 
sensors [17]. In this study, the reported computational savings were in the range of 59-77 percent. Full 
details of the automated substructuring algorithm are beyond the scope of this paper, and are briefly 
explained in this Section. In the developed algorithm, nonlinear criteria are set in each physics 
environment to identify the nonlinear regions in the FE model. These regions are isolated, and spectral 
techniques are used to automatically identify the locations of the linear regions in the FE mesh. The 
linear regions are then condensed in to a number of substructures, which are used later in the FE 
analysis. The generated substructures aspect ratio (AR) should be kept within a certain limit (1-10), in 
order to achieve accurate simulation results [18]. Therefore, the AR of the generated substructures is 
controlled within the specified limits using a modified version of the recursive spectral bisection (RSB) 
algorithm [19]. The stiffness of the generated substructures is fixed during the solution phase of the 
coupled field model. Therefore, only the stiffness of the nonlinear regions is updated during FE 
iterations, which considerably reduces the computational cost of the coupled FE analysis. Using the 
developed substructuring algorithm in the current study, computational savings ranging from 49-56 
percent were achieved. 

4. Validation Examples 

To assess the accuracy of the substructured multi-field coupled model, three examples were 
performed in this study. The first two examples are designed to validate the substructured EM-thermal 
model. The third example is utilized to validate the substructured structural-electrostatic model. 

4.1. Example 1 

In this example, the temperature distribution along the length of a fixed- fixed gold beam (4 = 400 
jim, Wm = 50 jum, tm = 2 jum), with an input power of 1 W is calculated. The beam was analyzed at two 
operational frequencies, 40 MHz and 40 GHz. At low frequencies, the current is spread evenly through 
the beam. For higher frequencies, self inductance of the conductor causes electrons crowding at the 
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outside edges of the beam, a phenomenon commonly known as the skin effect [5]. This phenomenon 
increases the heat dissipation to the beam; therefore higher temperatures are generated in the beam 
with higher operational frequencies. Results from the substructured coupled EM-thermal model 
highlight this effect, and are compared with numerical model data reported by Jensen et al. [5], in 
Figure 4. 

Figure 4. Substructured EM-thermal model predictions versus numerical model results 
reported by Jensen et al [5] of the temperature distribution versus coordinate along beam 
length for fixed-fixed gold beam at frequencies of 40 MHz and 40 GHz (/^ = 400 //m, 
Wm = 50 //m, tm^2 jum). 
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In their study, Jensen et al [5] modeled the switch in the upstate position using an iterative coupled 
field analysis. In their approach, the current density on the cross section of the beam is calculated using 
the finite element-boundary integral method. The calculated current density is then used to solve the 
steady state heat equation using 2-D finite element simulation. The two physics were calculated in an 
iterative manner until the solution converges to a constant temperature profile within 1% between 
iterations [5]. For more details, interested readers can refer to the study reported by Jensen et al. [5]. A 
good agreement between the two models is evident, as shown in Figure 4. 

4.2. Example 2 

In this example, the spatial average temperature rise versus input power in a fixed-fixed gold beam 
{Im = 400 jum, Wm = 20 jum, tm = 2 jum), is calculated. The beam was analyzed at operational 
frequencies of 2, 13.5 and 18 GHz. The results from the substructured EM-thermal model were 
compared with numerical and experimental data reported by Wang et al. [20]. The comparison is 
shown in Figure 5, where, EFE-BI refers to the finite element boundary integral method used to 
analyze the same problem by Wang et al. in 2006. 
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Figure 5. Substructured EM-thermal model predictions versus numerical model results and 
experimental data reported by Wang et al. [20] of spatially average temperature rise versus 
input power for fixed-fixed gold beam at frequencies of 2 GHz, 13.5 GHz and 18 GHz 

Qm = 400 jum, Wm = 20 jum, tm = 2 jum). 



T -T 




P(W) 



The results from the Major Vector algorithm were calculated initially based on a linearity 
temperature of 1 15 ^C. For an operational frequency of 13.5 GHz, the results start to diverge from the 
experimental measurements at an input power higher than 0.5 W. Therefore, the linearity temperature 
criteria were reduced to 95 for input power higher than 0.5 W, and the same criteria was used to 
calculate the spatial average temperature rise for frequencies above 13.5 GHz. As shown in Figure 5, 
the reduction of 7^, yields better solution accuracy, but this was on the expense of higher 
computational costs. 

4.3. Example 3 



In this example, results from the substructured structural-electrostatic model are compared with the 
closed- form analytical expression for the pull-in voltage given by [12]: 



(10) 



where Keffi^ the effective beam stiffness, and is given by Equation 11 for a uniform electrode [12]; 
is the permittivity of free space; ^4^^ is the effective area, given by Equation 12 [12,15,16], which 
mainly accounts for the fringing field effects and the charge redistribution effects: 



32E W f S.32N 

m mm _|_ ' ^ 



N = {a +(7^,)Wt 

V rp.s NL m \ 



(11) 
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m,ejj ^ m 



1 + 0.65 



W 



z 

P max 

So 



where, Gres = Ores{^-^m) IS the residual film stress; gnl= nEmZmaxl^lm^ '^^ an estimate of the induced 
axial stress due to the nonlinear stiffening; Zmax is the center deflection of the beam at a given applied 
voltage. 

The analytical and the FE simulation results, for the pull-in voltage versus beam length, for a fixed- 
fixed gold beam, at different microfabrication residual stresses are shown in Figure 6. The comparison 
shows a good agreement between the substructured structural-electrostatic model and the 
analytical solution. 

5. Results and Discussion 



(12) 



Given the good accuracy exhibited by the coupled multi-field FE model, the reliability of aluminum 
RF MEMS switches is examined in this section. Two candidate mechanisms are studied in this paper, 
mainly stiction problem at higher pull-in voltages and failure due to buckling. First, the effect of the 
operational frequency in the range of 0.1-100 GHz is studied, for flat switch membranes. Next, the 
effect of mean and gradient residual stresses, generated during different stages of microfabrication is 
evaluated. Finally, corrugated switch membranes and through holes in the switch membrane are 
presented, as alternative mechanical solutions for the RF switch reliability problem. The variation of 
the pull-in voltage {Vpi) and critical buckling strain {sc?) with switch length (/^) at different membrane 
thicknesses {tm = 0.5-3 //m), are shown in Figures 7 and 8, respectively. 

Figure 6. Substructured structural-electrostatic model predictions versus analytical results 
reported by Pamidighantam et al. [12] of pull-in- voltage versus beam length for fixed- 
fixed gold beam at different values of microfabrication residual stresses (fF^ = 50 jim, 
tjn= I jum, go = 2.5 jum , ^pi = 0.4, Em = 77 GPa). 
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Figure 7. Pull-in- voltage versus membrane length at different membrane thickness of flat 
aluminum RF MEMS switch {gm = = 0 MPa, Wm = 20 jum, T(x) = Tr), 
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Figure 8. Critical buckling strain versus membrane length at different membrane thickness 
of flat aluminum RF MEMS switch. 
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The pull-in voltage was calculated for a flat aluminum membrane, at room temperature, in the 
absence of residual stresses. It can be seen that Vpi decreases with the increase in 4 due to the 
reduction in the bending stiffness of the switch membrane [12]. The same effect can also be noticed for 
£cr, which varies proportionally with (tm/lm) . 

5.7. Effects of Operational Frequencies 



The steady state temperature distribution along switch membrane length at an input power of 
0.75 W, is plotted in Figure 9, as a function of operational frequency. 
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Figure 9. Temperature distribution along membrane length of flat aluminum RF MEMS 
switch at an operating frequency range of 0.1-100 GHz (Jm = 400 jum, Wm = 20 jum, 
t„ = 3^m,P = 0.75 W) 
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The maximum temperature occurs at the mid span of the beam. However, the maximum steady state 
temperature at 0.1 GHz (29.5 °C) is considerably lower than at 100 GHz (320 °C), due to the skin 
effect. The temperature rise induces compressive thermal stresses in the switch membrane, which 
increases the initial gap between the membrane and the bottom electrode (see. Figure 3). As a 
consequence, a coupled structural-electrostatic analysis of the switch membrane shows that Vpi 
increases dramatically with the increase in co, because of the change in the initial deformation state of 
the switch, as shown in Figure 10. Moreover, thermally induced compressive stresses cause the switch 
to buckle at CO > 10 GHz. 

Figure 10. Pull-in voltage versus membrane length of flat aluminum RF MEMS switch at 
operating frequencies of 10, 30 and 100 GHz (Wm = 20 jum, tm = 2 and 3 jum). 
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5.2. Effects of Residual Stresses 

Residual stresses are inherently generated during various stages of microfabrication. A general 
uniaxial residual stress field can be approximated by [21,22]. 

J? <Z<J? +t ^^^^ 

where, ctm, cfg/z) refer to the mean and gradient (intrinsic) residual stresses generated in the thin film 
during microfabrication, respectively. Typical values reported in the literature indicate experimental 
measurements in the range of 0 to over 150 MPa for au and 5-30 MPa for aGr(z) [8]. Uniaxial tensile 
mean residual stresses (50-100 MPa), were applied to the flat aluminum RF switch to examine its 
effect on Vpj and average axial strain (s). The pull-in voltage-frequency curves predicted by the multi- 
field coupled model are shown in Figure 11, at different values of om for flat RF switch membrane 
with dimensions (Im = 600 jum, Wm = 20 jum, tm = 3 jum). 

It can be seen that the applied residual stresses increase the pull-in voltage by 20^0 percent. This 
increase can be explained by the increased stiffening of the beam due to the increase in the value of the 
applied residual stress. However, as shown in Figure 12, the tensile residual stresses counteract the 
compressive thermal stresses induced at higher operational frequencies. 

It can be seen that s can be shifted above the value of the critical buckling strain, which increases 
the power handling capability of the switch, to higher operational frequencies. Therefore, the 
microfabrication tensile residual stresses can be utilized to eliminate the buckling failure. Conversely 
the stiction problem, due to charge build-up in the dielectric layer associated with higher actuation 
voltages, can still cause switch failure at a limited number of operational cycles. 

Figure 11. Effect of microfabrication mean residual stresses on Vpi of flat aluminum RF 
MEMS switch at operating frequency range of 10-100 GHz (4 = 600 jum, Wm = 20 jum, 
tm = ^ jum ). 
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Figure 12. Effect of microfabrication mean residual stresses on s of flat aluminum RF 
MEMS switch at operating frequency range of 10-100 GHz (4 = 600 jam, Wm = 20 //m, 
= 3 iim). 
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5,3. Mechanical Approaches 



As discussed in subsection 5.2, a proper control of tensile residual stresses generated during 
microfabrication, can be utilized to increase the power handling capabilities of RF MEMS switches. 
However, the higher actuation voltages exhibited; because of membrane stiffening is still a major 
problem. Therefore, to achieve reliable switch designs, much lower actuation voltages should be 
achieved. Three strategies have been proposed in the literature to achieve this goal. The first one 
consists in modifying the out-of-plane geometry of the switch membrane [13]. The second strategy 
consists in the use of corrugated switch membranes [13]. The third strategy consists in the addition of 
holes to reduce the membrane effective Young's modulus [22]. The second and third strategies are 
more achievable in microfabrication [13]; hence, the combination of both strategies in RF MEMS 
switch design is examined in this subsection, for an RF MEMS switch with dimensions (Im = 600 jum, 
Wm = 20 jum, tm = 3 jum). 

First, the effect of number of corrugations (Nc) and corrugation thickness (tc), expressed by 
corrugation thickness ratio (tc/tm), was investigated. The pull-in voltage versus Nc is shown in 
Figure 13. This analysis was conducted for ctm range from 50-100 MPa, tc/tm = 0.16 and co = 30 GHz. 
It can be seen that a 20 percent reduction in Vpi was achieved with Nc = 24, 

The pull-in voltage versus tc/tm is shown in Figure 14. This analysis was conducted for ctm range 
from 50-100 MPa, Nc =24 and co = 30 GHz. It can be seen that proper selection of tc can dramatically 
reduce Vn, up to 60 percent, at the studied range of microfabrication residual stresses. 
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Figure 13. Effect of number of corrugations on Vpi of corrugated aluminum RF MEMS 
switch at ctm range of 50-100 MPa (/^ = 600 fim, Wm = 20 jum, tm = 3 jum, tc/tm = 0.16, 
CO = 30 GHz). 
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Figure 14. Effect of corrugation thickness ratio on Vn of corrugated aluminum RF MEMS 
switch at om range of 50-100 MPa (Im = 600 jum, Wm = 20 jum, tm ^ ^ jum, Nc = 24, 
CO = 30 GHz). 
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The average axial strain was also examined for RF switch (Nc = 24, tc/tm = 0.5), at range from 
10-100 GHz, as shown in Figure 15. It can be seen that for aM= 50 MPa, the buckling failure occurs at 
CO ^ 50 GHz, while for aM= 100 MPa, the power handling capability of the RF switch can be increased 
beyond 100 GHz. 
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Figure 15. Average axial strain versus cw of corrugated aluminum RF MEMS switch at om 
range of 50-100 MPa {Nc = 24, tc/tm = 0.5). 
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A very promising design is the addition of holes in to corrugated switch membranes [4,22]. The 
holes total area can be up to 60 percent of the total surface area of the switch structure [4]. The 
perforation pattern is mainly characterized by the ligament efficiency, //// = In/pitch, defined as the 
ratio of the remaining link width to pattern pitch, as seen in Figure 16. 

Figure 16. Parameters for the ligament efficiency in RF MEMS switch membrane with 
through holes. 
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The holes allow for the release of part of the residual stresses in the membrane, and reduce the 
effective Young's modulus of the switch membrane. The reduction of the effective Young's modulus, 
yields lower values of actuation voltages. Moreover, this reduction increases the positive residual axial 
strain component, which counteracts the compressive strains generated at higher frequencies due to the 
skin effect. Therefore, the power handling capability of the RF switch can be increased considerably, 
without buckling, while still at reliable ranges of actuation voltages. Figure 17 shows the average axial 
strain versus an operational frequency range 0-200 GHz. 
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Figure 17. Effect of ligament ratio (jih) on average axial strain of corrugated aluminum RF 
MEMS switch at co range of 30-200 GHz {om = 100 MPa, Nc = 24, tc/tm = 0.5) 
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This analysis was conducted for juh range from 0.4-0.8, gm = 100 MPa, Nc = 24 and tc/tm = 0.5. It 
can be seen that a higher operational frequency range is achievable for corrugated switch designs with 
holes (co > 200 GHz, juh = 0.4), compared to ~ 100 GHz, for corrugated switch design. 

6. Conclusions 

The commercialization of RF MEMS switches is subject to the elimination of the buckling failure at 
high RF power and the stiction problem associated with higher actuation voltages. In this paper, a 
substructured coupled multi-field model was presented to examine the performance of RF MEMS 
switches at an operational frequency range 0.1-100 GHz. The effect of microfabrication residual 
stresses and effectiveness of some mechanical approaches has been examined in detail. The coupled 
model, proposed in this study, represents the first effort to link an actual temperature distribution, due 
to the skin effect, associated with higher operational frequencies to a coupled structural-electrostatic 
analysis of RF MEMS switches. Other studies reported in the literature assume a uniform temperature 
distribution in the switch membrane, which might generate inaccurate predictions of the actual pull-in 
voltages. The main conclusions from the analyses conducted in this study may be summarized 
as follows: 

1. The temperature rise at higher operational frequencies induces compressive stresses in the 
switch membrane of flat RF MEMS switches. The induced compressive stresses lead to buckling 
and device failure for co > 10 GHz. Moreover, the change in the deformation state of the switch 
leads to a dramatic increase in the required actuation voltages. 

2. In-plane tensile residual stresses generated during microfabrication counteract the induced 
compressive thermal stresses. In this sense, proper control of microfabrication residual stresses can 
increase the power carrying capacity of the RF switch. However, higher values of actuation voltages 
are exhibited because of the membrane stiffening effect. 
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3. Membrane corrugations, at a distance from the support area can be used to maintain the 
actuation voltages at lower values. Therefore, using these membranes in conjunction with a proper 
control of microfabrication residual stresses, increase the reliable range of operation of the RF 
switch. 

4. A better design can be achieved by introducing holes to corrugated switch membranes. These 
holes help in reducing the effective Young's modulus of the switch membrane. This mainly would 
yield lower actuation voltages. Moreover, this reduction increases the positive residual axial strain 
component, which counteracts the compressive strains generated at higher frequencies. Therefore, 
the power handling capability of the RF switch can be increased considerably, without buckling, 
while maintaining a suitable actuation voltage. Our analysis of a corrugated switch membrane, with 
holes at a ligament efficiency of 0.4, shows that this design can increase the reliable operation of the 
RF switch to (D > 200 GHz. 
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